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Abstract 

It is shown that, as functions of their rank, cumulants of QCD 
multipHcity distributions in small phase space bins possess the quasi- 
oscillating behavior similar to that found for them in the total rapidity 
range. First minimum moves to lower ranks for smaller bins. For the 
total rapidity range, it moves to higher ranks with energy increase. 
The running property of the QCD coupling constant is in charge of 
these effects which can be verified in experiment. 

Recently, it was shown [1] (for a review, see [2]) that cumulants of multi- 
plicity distributions in QCD behave in a rather peculiar way revealing pro- 
nounced oscillations. The analysis was done for the total rapidity range. 
Its predictions have been confirmed by experiment [3]. From the theoretical 
side, those results imply that QCD distributions do not belong to the class of 
infinitely-divisible ones (in particular, they can not be fitted by the negative 
binomial distribution - NBD) and prohibit Poissonian cluster models [4]. 
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Here, we shall answer two more questions related to the behavior of cu- 
mulants. First, we show that the first minimum of the ratio of cumulant to 
factorial moments for the total rapidity range moves to higher ranks when 
the primary energy increases. Second, we show that the similar minimum ap- 
pears when multiplicity distributions in small phase space bins are analyzed^ 
and it shifts to lower ranks for ever smaller bins. 

The above properties can be easily guessed from the formulae obtained 
in the original paper [1]. Therefore, by using them, we shall, first, give quali- 
tative arguments in favor of such properties, then show how different factors 
influence the conclusion, and, finally, discuss its experimental implications. 

In its simplest double-logarithmic approximation (DLA), quantum chro- 
modynamics fails to describe experimentally observed multiplicity distribu- 
tions predicting very wide shape of them. However, already in the next 
(modified leading logarithm - MLLA) approximation its predictions become 
quite reasonable. Besides considering just the multiplicity distributions P„, it 
is convenient to use their factorial (Fq) and cumulant (Kg) moments defined 
by the relations 

_ {n{n - l)...{n - g + 1)) _ Y.njn- l)...{n - g + 1)P„ 

9-1 

Fq = ^ C,^ Kq-rnFm-i (2) 
m=0 

where are the binomial coefficients, and the set of recurrence relations 
(2) is applied to define the cumulants when the factorial moments are deter- 
mined according to eq.(l). 

As a function of the rank g, their ratio 

H, = ^ (3) 

behaves like in DLAQ but acquires a minimum in MLLA crossing once 
the abscissa axis. Higher-order approximations predict further crossings and 

^Such an analysis has become very popular in connection with ideas of intermittency 
and fractality (for a review, see [5]). 

•^To get an insight, let us notice that the negative binomial distribution with a param- 
eter k interpreted as a number of independent sources predicts the ratio H^^^ = q^^ . 
Herefrom, one concludes that DLA corresponds to NBD with rather small value of the 
parameter k = 2 and, consequently, gives rise to very wide multiplicity distribution. 
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quasi-oscillatory behavior of the ratio (see [2]). The position of the first 
minimum, already observed in MLLA, seems rather stable. It depends on 
the QCD coupling constant only and is not strongly influenced by the higher- 
order terms. In gluodynamics, it is given [1] approximately as 

?mm~ 7:^ + ^ + 0(70), (4) 

where 70 = (Gas/vr)^/^, is the QCD coupling constant, and hi = 11/24 is 
the coefficient calculable in MLLA. Higher approximations contribute terms 
0(70), and, in practice, can be neglected when dealing with the position of 
the very first minimum. 

The formula (4) immediately shows to us that the position of the first 
minimum moves to higher values for jets with larger virtuality Q since the 
QCD coupling constant is running as ln~^ Q. Therefore, the predicted shift 
of the minimum is 

Qmin OC lu^^^ Q (5) 

with increase of virtuality Q. Thus, analyzing the multiplicity distributions 
at ever higher energies one would notice the shift of the minimum of Hq to 
higher rank values. 

One could guess already that the position of the minimum would shift to 
smaller values of qmin if multiplicity distributions within small phase space 
bins are considered. Really, it would mean considering the sub jets with lower 
virtuality, and, therefore, with larger effective values of as and 70. To confirm 
this statement and estimate the various factors, contributing to it, one should 
perform the analysis anagolous to that done in [6] where factorial moments 
for small phase space bins were treated. Referring for further details to [6], 
let us outline the main steps of the analysis. 

To calculate the moments of the distributions within small bins, one 
should convolute three stages of the process of the emission of a gluon jet by 
a quark: 

1. the initial emission of a hard primary gluon jet by a quark, 

2. its subsequent evolution giving rise to secondary jets with some energy 
spectrum. 
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3. the choice (among those jets) of the final subjet which hits the bin 
under consideration and the estimate of its correlators. 



It has been done in [6] for factorial moments of the distributions in small 
phase space bins but can be repeated with slight modifications for cumulants 
as well. For example, for cumulants one writes down the formula analogous 
to (3.2) of [6] : 



e y ' " J E 2n ^ \QJ J k \k 

{n{ke;l))^H,{kQ-l), (6) 

where the energy E integration with Altarelli-Parisi kernel $ corresponds to 
the stage 1), the energy k integration with the spectrum is according 
to the stage 2), and the factor {n)'^Hg describes the cumulant of the subjet 
hitting the angular window 6 (the stage 3)). Here Q = EqQq,Eq is the 
primary jet energy and Bq ~ 1. The only difference from the integrals for 
factorial moments in [6] is due to the factor Hg for the final subjet which 
should be evaluated. However, this factor is a slowly varying function of 
k while the spectrum and the g-th power of multiplicity change very 
rapidly. Therefore, the steepest-descent method can be again applied to cal- 
culate the correlator (6) with Hg estimated at the stationary point. Since the 
unnormalized factorial moments are given by the similar expressions without 
Hg in the integrals, one gets, after dividing cumulants by factorial moments, 
that the ratio Hg for the multiplicity distribution in a small bin is just the 
ratio for the subjet hitting the bin when calculated at the stationary point 

= {kQ)s.p.- Because the subjet products cover completely the bin 9, its 
ratio Hg is to be calculated from the same equations as for the total rapidity 
rangeQ and, therefore, possesses the similar quasi-oscillatory behavior. How- 
ever, the positions of maxima and minima are shifted because it is calculated 
now at the stationary point but not for the total range of angles. 

We shall deal with the behavior of the first minimum which is well de- 
scribed by MLLA where the ratio Hg is given [1] by the equation 

rjMLLA 7o(l-2/lig7) fr,. 
= ^2 • 



^The only difference is that the total phase space should be replaced by the phase space 
window under consideration. 
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Here, 7 is an anomalous dimension of QCD determining the energy depen- 
dence of mean multiplicity by the relation (n) = exp^J dtj{t)). The only 
problem now is to insert 7 for small bins and to use the stationary point 
value of 7o. Actually, the value of 7 in MLLA for small bins was calculated 
in [6] (see eq.(3.30) there) and, in case of gluodynamics, is given by 

7 = 7o(l-Cq7o), (8) 

' 48 + + • ^ ' 

One easily gets that in the region of interest to us g = 3 — 5 the correction 
term is negligible 

Cg7o < 1- (10) 

Therefore, for our purposes 7 ~ 70, and one gets the result (4) with high 
enough precision. Thus, MLLA modification of the anomalous dimension 7 
is not important here. 

The stationary point of the integrand in eq.(6) is estimated using the stan- 
dard methods [6, 7]. For qualitative estimates, we are aimed to, it is enough 
to demonstrate the behavior of 70 at the stationary point as calculated for 
fixed coupling since influence of its running on the position of the stationary 
point is a higher-order effect. Then it is easy to get for the virtuality Kq of 
the sub jet at that point 

.„ = «(g^) . (11) 

i.e. "the effective mass" of the subjet diminishes almost proportionally to 
the opening angle with additional (and not important for our purposes) 
factor 0V9' appearing due to the energy spectrum D® decrease. Therefore, 
the slight shift of the minimum position q„iin to smaller values is envisaged 
at small bins since g„jj„ oc contains the term ~ In^^^ Kq (let us note that 
kq replaces Q compared to eq.(5)). 

Even though more rigorous expressions are available, we decided to use 
the semiquantitative description due to two reasons. First, mathematical 
details are rather cumbersome but do not change main conclusions presented 
above. Second, the experimental verification of the results is also possible 
at the semiquantitative level only so that the theoretical precision should 
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not exceed the experimental facilities. That statement is easy to appreciate 
if one reminds that experimentalists can get the ratio Hq at integer values 
of q only. They proceed from the factorial moments evaluated according to 
eq.(l) to the cumulants defined through the set of recurrence relations (2) 
for integer points q . Therefore, they can not observe the continuous drift of 
qmin with energy or with bin sizes. Nevertheless, such a shift (if noticed as 
a shift of the minimum position from one integer rank value to another) will 
be of interest since it demonstrates further details predicted by QCD. 

Probably, the best way to find it out in experimental data is to study 
near the minimum point qmin the evolution of the differences Fq — Fq^i with 
increase of energy and/or with decrease of bin sizes. Unfortunately, it is hard 
to measure high-rank cumulants at small bins with high enough precision. 
Besides, they are very sensitive to the selection procedure in a given exper- 
iment (even though this sensitivity can be, probably, used to compare the 
criteria attempted by different experimental groups). 

In conclusion, we have shown that the minimum of the ratio of cumulant 
to factorial moments of multiplicity distributions predicted by QCD [1] and 
confirmed by experiment [3] is subject to the shift to higher values of ranks at 
higher energies and to smaller ranks when multiplicity distributions in small 
phase space bins are measured. Experimental verification of those effects 
would demonstrate once again the running property of the QCD coupling, 
and, furthermore, the very applicability of QCD to description of subtle 
correlation effects even for rather soft processes of multiparticle production. 
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